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Abstract. We present the hard X-ray time-averaged spectrum of the intermediate polar V709 Cas observed with 
INTEGRAL. We performed the observation using data from the IBIS/ISGRI instrument in the 20-100 keV energy 
band and from JEM-X at lower energy (5-20 keV). Using different multi-temperature and density X-ray post- 
shock models we measured a improved post-shock temperature of ~ 40 keV and estimated the V709 Cas mass to 
be 0.821q;25 Mq. We compare the resulting spectral parameters with previously reported BeppoSAX and RXTE 
observations. 
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1. Introduction 

V709 Cas (RX J0028+5917) is a member of the inter- 
mediate polars (IP) systems, a sub-class of magnetic cat- 
aclysmic variables. This system consists of an accreting 
white dwarf (WD) and a low-mass late type main se- 
quence companion star. The accretion onto V709 Cas 
is believed to be driven through the Roche-lobe over- 
flow, where the accretion flow from the secondary pro- 
ceeds towards the IP through an accretion disk, until it 
reaches the magnetospheric radius. Here the material at- 
taches to the magnetic field lines and follows them al- 
most radially at free-fall velocity towards the magnetic 
poles of the IP surface. At some distance from this sur- 
face, the accretion flow undergoes a strong shock, below 
which material settles onto the IP, releasing X-ray as it 
cools by thermal bremsstrahlung processes (see reviews 
lAizul Il973t ICropperl Il99fi IPattersonl Il994|) . According 
to this standard model the spectra of the X-ray post- 
shock emitting reg ion has a multi-temp erature and den- 
sity structure (e.g.. lCropper et aTl ll999'). In order to mea- 
sure the maximal shock temperature, hight energy obser- 
vatories like INTEGRAL are needed. In the most case 
the temperature of the post-shock region of the accre- 
tion column of IPs arc in the order of ~ 10-60 keV. A 
hard X-ray study can also be used to estimate the IP 



mass by measur ing the maximum temperature of t h e post - 
shock plasma i Suleimanov. Revnivtsev fe Ritterl l2005|) . 
Using the INTEGRAL / RXTE data on the IP V1223 Sgr 
iRevnivtsev et al.l l)2004|) studied such a broadband X-ray 
spectrum to determine the shock parameters and the IP 
mass. 

The X-ray V709 Cas source was discovered and 
identified as a IP fro m the ROSAT All sky Survey 
l|Haberl fc Motchl. 119951) and has been extensively s t udied 
fTaherl fc Motchl 1199.4 iMotch et all [1 996; 



with ROSAT 
iNorton et all 



2|X 

, 199 9J), with a joint RXTE /BeppoSAX X 

ray observation (|De Martino et all l200ll) and using opti 



cal spectroscopy and phot ometry (|Bonnct-Bidaud e t all 
l200lt iKozhevnikwl l200l|) . resp ectively. The pu ls e pe- 
riod was found to be 312.8 s l)Haberl fc Motchl. Il995t 
INorton et all Il999j) and a detailed spectroscopic study 
finally concluded to a 5.34 h orbital period after 
a previous ambiguity bet ween 5.4 hr and 4.5 hr 
l|Bormet-Bida,ud et al 1 Eiol . 

We present here a first INTEGRAL results of a broad- 
band spectrum study, from 5-100 keV, where the energies 
above 20 keV are based on hard X-ray IBIS/ISGRI obser- 
vations of V709 Cas. 



2. Observations and Data Analysis 
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The present dataset was obtained duri ng the Target 
of Opportunity (ToO) A02 INTEGRAL l|Winkler et all 
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Fig. 1. The 20-40 keV IBIS/ISGRI mosaicked and de- 
convolved sky image of the ~ 182 ks observation. Image 
size is ~ 12° x 5°, centered at V709 Cas position. The pixel 
size is 5'. V709 Cas and IGR J00291+5934 at a separation 
of 20' were detected at a significance of ~ 12.6<7 and at 
~ 82.8cr, respectively. 



2003) observation of the Cas A region performed from 5- 
6 and 7-9 December 2004 (53345.6-53346.8 and 53347.8- 
53349.8 MJD), i.e. from part of INTEGRAL satellite rev- 
olutions 262 and 263. We use data from the IBIS/ I SGRI 
(20-100 keV) coded mask imager (|Ubertini et all 120031: 
iLebrun et al. , l2003h for a total expos ure of 181.9 ks and 
from the JEM-X (5-20 keV) monitor l|Lund et all [2003) 
for a total exposure time of 82 ks. For ISGRI, the data 
were extracted for all pointings with a source position 
offset < 7°, and for JEM-X with an off set < 3.5°. The 
spectrometer (SPI) l|Roaues et all |2003) was not used to 
extract the hard X-ray spectrum due to the lower sensi- 
tivity of this instrument with respect to IBIS/ISGRI for a 
weak source below 100 keV. Above 90 keV, V709 Cas was 
not consistently detected in a single exposure and in the 
total exposure time. Data reduction was performed using 
the standard Offline Sci ence Analysis (OSA) version 4.2 
l|Courvoisier et aLll2003l) . The algorithm s used in the spa- 
tial an d spectral analysis are described in lGoldwurm et all 
(20031). 



3. Results 

3.1. ISGRI Imaging 

Fig. ^ shows a significance map around the source 
V709 Cas in the 20-40 keV energy range. Single pointings 
were deconvolved and analyzed separately, and then com- 
bined in mosaic images. Two sources are clearly detected 
at a significance level of respectively 12.6(7 for V709 Cas 
and 82.8cr for IGR J00291+ 5934, a newly discove red mil- 
lisecond pulsar in outburst ()Falanga et alll2005|) . In the 
energy band 40-80 keV, the confidence level was 5(7 for 
V709 Cas and 46cr for IGR J00291+5934. At higher ener- 
gies above 100 keV V709 Cas was not detected at a statis- 
tically significant level cither in single exposures or in the 
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Fig. 2. The unfolded combined INTEGRAL JEM- 
X/ISGRI (5-100 keV) spectrum of V709 Cas with the 
best- fit post-shock model (line). Residuals between the 
data and model are shown in the bottom panel in units of 
sigma. 



total exposure time. To obtain precise source locations we 
simultaneously fitted the ISGRI point spread function to 
the two close sources. We obtained a position for V709 Cas 
at ajaooo = 00 h 28 m 55!29 and feooo = 59°16T4'.'0. 
The position of IGR J00291+5934 is given by a J20 oo = 
00 h 29 m 02!92 and <5 J2000 = 59°34'06'.'4. The source posi- 
tion offsets with respect to the optical catalog positions 



llDownes. Webbink fc Sharalll997l:lFox fc KulkarnAl2004h 
are 1.'4 for V709 Cas and 0.'2 for IGR J00291+5934. The 
errors are I'.b and 0.'2 for V709 Cas and IGR J00291+5934, 
respectively. These are within the 90% confi dence level 
assum ing the source location error given by lOros et alJ 
(2003J). The derived angular distance between the two 
sources is ~ 20'. Due to the fact that INTEGRAL is 
able to image the sky at high angular resolution (12' for 
ISGRI and 3' for JEM-X), we were able to clearly dis- 
tinguish and isolate the high-energy fluxes from the two 
sources separately. This allows us to isolate and study 
the X-ray emission of V709 Cas during the outburst of 
the IGR J00291+5934 pulsar. The lack of contamination 
by the pulsar outburst was verified by building V709 Cas 
light curves in different energy ranges. While during the 
observation, the pulsar showed a significant decay with 
an e-folding time of ~ 6.6 day ()Falanga et alll2005l) . the 
V709 Cas (20-80 keV) light curve remains stable with a 
constant mean counting rate of ~ 0.6 ± 0.09. 

3.2. Spectral Analysis 

The spec t ral an alysis was done using XSPEC version 11.3 
l|Arnaudl Il996|) . combining the 20-100 keV ISGRI data 
with the simultaneous 5-20 keV JEM-X data. Due to the 
short exposure time of JEM-X, and therefore lower statis- 
tics, we rebinned this data in a 5 channel energy response 
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matrix. A constant factor was included in the fit to take 
into account the uncertainty in the cross-calibration of the 
instruments. A systematic error of 2% was applied to the 
JEM-X/ISGRI spectra which corresponds to the current 
uncertainty in the response matrix. All spectral uncertain- 
ties in the results are given at a 90% confidence level for 
single parameters. 

The joint JEM-X/ISGRI (5-100 keV) broadband spec- 
trum was first fitted with a simple optically thin thermal 
bremsstrahlung. The data fitted with x 2 /dof=5/12, and 
the plasma temperature was found around 27 keV. In or- 
der t o compare with previou sly reported measurements 
fe.g.. lDe Martino et allEoOlj) . we also fit the plasma emis- 
sion model Mekal. We found a consistent plasma temper- 
ature value with the previous thermal bremsstrahlung fit, 
where the flux in the 5-100 keV energy range is 1.2 xlO -10 
erg cm~ 2 s _1 (see TableQ}. During our INTEGRAL obser- 
vation the source flux in the 3-100 keV was ~ 1.4 times 
smaller than that observed in July 1 998 by BeyyoSAX 
and ~ 2 times smaller by RXTE foe Martino et all 
l200li ISuleimanov. Revnivtsev fc Ritterl l2005[h The ob- 
tained value of the temperature parame ter (26 keV) is sig- 
nifican tly smaller than that obtained bv lDe Martino et all 
l|200l[) from BeppoSAX data ( 42 keV) and RXTE data , 
(36 keV) using the same model, be Martino et all l|200l[) 
found the same plasma temperature ~ 26 keV, only by in- 
cluding a reflection model in the fit. We attempted to test 
the hypothesis of a Compton reflec tion in the spectrum in- 
cludin g the reflection model from iMagdziarz fc Zdziarski 
l)l995|) . The lack of statistics in the low energy part of the 
JEM-X data prevent us to determine significant reflection 
parameters. The best fit parameters using only the hight 
energy ISGRI data alone agree with the value from the 
joint JEM-X/ISGRI spectrum. This confirm the impor- 
tance of the hight energy observation to determine the 
post-shock temperature above 20 keV. 

However, from the standard accretion shock model, 
when the matter falls onto the WD, the X-ray post-shock 
emission region is expected to show a multi-temperature 
and multi-density structure. To fit the spectrum we there- 
fore use a more physically motivated model where the X- 
ray emission is determined through a density and temper- 
ature gradient along the emission region. The broadband 
time-averaged spectrum is then given by summing local 
bremsstrahlung spectra in the region between the WD sur- 
face and the shock dista nce, zp. We use the same geometri- 
cal mo del described bv ISuleimanov. Revnivtsev fc Ritter 
where the observed flux is given bv IZombeck 



119991) . In the flux equation, /i = 0.62 is the mean molecular 
weight of fully ionized accreting matter, toh is the hydro- 
gen mass and k is the Boltzmann constant. The only two 
parameters in this model is now the WD mass, M w( i, and 
the accretion rate by unit of surface a. From the best fit, 
using a standard local accretion rate of a = 1.0 g cm -2 
s _1 , we found M w d = 0.82 M Q corresponding to a shock 
temperature of 39 keV (Tabled . The % 2 value is compara- 
ble with one of the simple bremsstrahlung model, however 
the post-shock model represent in more detail the physical 
processes in the WD emission region. Absorption by neu- 
tral hydrogen and partial absorber have not been included 
in these models since they have no significant effect above 
5 keV. TableQ]gives the best fit parameters of this column 
model, for the ISGRI and combined JEM-X/ISGRI data 
sets. In Fig. [3 we present the vF v spectrum of the entire 
observation, plotted together with the residuals in units 
of a with respect to the best fit post-shock model. 

3.2.1. Spectral Analysis: Models 

Different models for the structure of accretion column 
have bee n put forward based on assumptions of constant 
press ure IIFrank et al. , 2002), influence of the magnetic 
field ilWu. Chanmugam fc Shavhl Il994) or gravitational 
potential ijGropper et al.lll999h . In an effort to also eval- 
uate the influence of these different assumptions, a phc- 
nomcnological fit with a variable emission measure can 
be performed. Each model is indeed characterized by a 
power-law type density and temperature profiles along the 
accretion column. Assuming profiles with (T/T s h oc k) = 
(x/x shock ) a and (p/p s hock) = (z/^shock)' 3 , it can be shown 
that the emission measure EM — p(z) 2 Adz is de- 
fined as: 



EM = {T/T shock y with r = 



_ (2/3+1) 



(2) 



1990): 



9.52 x 10~ 38 



T(z) 



-1/2 



E 



kT{z) 



-0.4 



exp 



E 



kT(z) 



dz, (1 



allowing the density profile, p(z), and temperature pro- 
file, T(z), to vary in the post- s hock r egion (e.g., see 
ISuleimanov. Revnivtsev fc Ritterl 1200.4 ICronper etail 



The EM therefore also follows a power-law in tempera- 
ture with an index T and the column spectrum can be de- 
scribed with a power-law multi-temperature plasma emis- 
sion model such as the XSPEC-CEMKL. The parameters 
a and /3 have been determined by fitting the published 
profiles of the different accretion shock structure models 
and are reported in tabled To determine the post-shock 
temperature, the JEM-X/ISGRI spectrum have been then 
fitted using CEMKL and fixing the model dependent power- 
law index, T, and the results are reported also in table 

4. Conclusions 

A first broadband spec trum result on V709 Cas 
was already studied by foe Martino et all l200lj) us- 
ing BeppoSAX and RXTE data. We report here the 
first 5-100 keV high energy INTEGRAL spectrum result 
of V709 Cas. We found a Bremsstrahlung temperature 
around 26 keV which is in agreement of a post-shock tem- 
perature around 40 keV (see Fig.[3J). The theoretical spec- 
trum of the post-shock model shown in Fig.^is calculated 
using equation (1), where the bremsstrahlung spectrum is 
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Table 2. Best fit parameters from the JEM-X/ISGRI data in the 5-100 keV energy range. 



Model 


a 




r 


B-Field (MG) 


fcT s hock (keV) 


X 2 /dof 


Fx (erg cm 2 s 1 ) 


Frank et al. 2002 
Wu et al. 1994 
Wu et al. 1994 
Suleimanov et al. 2005 


0.4 
0.341 
0.389 
0.312 


-0.4 
-0.399 
-0.452 
-0.433 


0.5 
0.592 
0.247 
0.430 


30 


39.021 JJ;2 

AO q— 13.8 

40.6-3 2 :? 


4.91/12 
4.8/12 
5.07/12 
4.95/12 


1.36x10" 10 
1.35x10- 10 
1.4xl0" 10 
1.37x10""' 



Table 1. Best fit parameters of the phase-averaged X-ray 
spectra. 



Dataset 
Energy range 


JEM-X/ISGRI ISGRI 
(5-100) keV (20-100) keV 


Model 


Bremsstrahlung 


kT Blems (keV) 
X 2 /dof 

F x (erg cm 2 s x ) 


2fi 7" 6 - 7 25 5+ 9 - 3 
5.3/12 5.8/9 
1.2X10 -10 0.52x10""' 


Model 


Mekal 


fcT M ax (keV) 
X 2 /dof 

F x (erg cm 2 s 1 ) 


25.7±l:l 24.6l 8 5 6 s 
5/12 5.8/9 
1.3X10 -10 0.57x10""' 


Model 


Post-Shock 


M wd (M ) 
^l ock (keV) 
X 2 /dof 

F x (erg cm 2 s 1 ) 


0.82«i 2 0.85±° ;« 
39.38 41.3 
5/12 6.3/9 
1.2xl0- 10 0.95x10-"' 



Calculated. 



the best fit using the same equation as a bremsstrahlung 
model. The Bremsstrahlung temperature is lower then 
Tshock since TB re ms represent more a weighted mean value 
between the temperature distribution from the IP surface 
to the shock hight z . Using JEM-X data in the energy 
band 5-20 keV and statistically significant ISGRI data 
in the energy range from 20-100 keV we estimated the 

M w . Using the broad band 



1+0.12 



source mass to by 0.82_ n ^ 

RXTE spectrumlSuIe imanov. Revnivtsev fc Ritterl l)2005|) 
found within the error bars the same IP mass. The maxi- 
mum shock temperature reported in tabled is calculated 
from: 

_ 3 fx G Mwd mH 

-L shock — / d 

O & Kind 

Using different accretion shock structure models we found 
the same post-shock temper ature (see ta ble HI). The IP ra- 
dius was calculated from the lNauenberd l|l972j) WD mass- 
radius relation and we obtained a radius of R w d = (0.68 ± 
0.13) x 10 9 cm. The obtained estim ation of the IP mass is 
slightl y larger than that obtained bv lBonnet-Bidaud et ID 
(2001). These authors inferred the mass estimation using 
optical spectroscopy, i.e. Balmer absorption lines, believed 
to be from the WD photosphere. 

We estimates the accretion rate from the X-ray lumi- 
nosities of the INTEGRAL observation. For the source dis- 



tance we use the range 210-250 pc l|Bon nct-Bida ud et all 
1200 111 . The absorption corrected broad band X-ray flux of 
V709 Cas in the 0.5-100 keV energy range is F . 5-100 koV = 
1.6 x 10" 10 erg cm" 2 s" 1 , therefore the isotropic luminos- 
ity of one accretion colum is L Xi i — 47rZ? 2 Fo. 5-100 keV = 
0.84— 1.2 x 10 33 erg s _1 . From the Bremsstrahlung model 
the emission measure of one emitting accretion colum is 
E mA = J Nl dV = 4.72 x 10 12 4irD 2 cm" 3 , where N e is 
the number density of the electron density in a hot emit- 
ting plasma. The volume emission measure therefore is 
E,, 



found bv lDe Martino et all (|200l|h As the X -ray emission 
is the major channel of energy losses in the post-shock 
regio we derived the mass accretion rate M 

M ~ 3 At m H G2L Xi i 



8 kT t 



hock 



The factor 2 stands in order to account for two emitting 
magnetic poles of the IP. Substituting measured values we 
obtained M = 0.82 - 1.2 x 10 16 g s" 1 . Using m wd = 0.82 
in units of the solar mass, R w d — 0.68 x 10 9 cm and the 
M we evaluated the height of the standing shock h s and 
the number density of elec trons in the post-shock region 
N e following IWarnerl l|l995lh 



N P = 3.1 x 10 15 M 



16 



h s = 9.6 x 10 7 m wd 



m wd 



N K 



1/2 



R 



-3/2 



10 16 cm- 3 



i? 9 1 cm, 



,-3 



where Mie is the mass accretion rate of 10 g s , Rg 
is the IP radius in the units of 10 9 cm, and / is the 
fraction of the IP surface, occupied by two accretion 
columns. We derived N e = 4.1 — 5.8 x I0 15 fl^ cm" 3 and 
h a = 0.21 - 0.28 /_ 3 R wd , where /_ 3 = 10" 3 / is the 
fraction of the surface of the IP, occupied by two accre- 
tion columns in units of 10~ 3 . The derived parameter s 
are similar to those obtained bv lDe Martino et al . (2001); 
iBonnet-Bidaud et alJ (|200l|) . 
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Fig. 3. The figure shows the agreement between the found 
value of the shock temperature from the post-shock model 
(see Sec. 3.2) with the best fit value using equation (1) as 
a Bremsstrahlung model. 
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